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In the framework of nonextensive statistical mechanics, the equilibrium structures of astrophysical 
self-gravitating systems are stellar polytropes, parameterized by the polytropic index n. By careful 
comparison to the structures of simulated dark-matter halos we find that the density profiles, as 
well as other fundamental properties, of stellar polytropes are inconsistent with simulations for any 
value of n. This result suggests the need to reconsider the applicability of nonextensive statistical 
mechanics (in its simplest form) to equilibrium self-gravitating systems. 

PACS numbers: 



I. INTRODUCTION 

Nonextensive statistical mechanics is a generalization 
of thermodynamics and statistical mechanics proposed 
by C. Tsallis in 1988 aiming to overcome the limita- 
tions of the Boltzmann entropy in its conventional appli- 
cations 0, [H, 0, H[ . The theory has considerably widened 
the fields of application of statistical physics, allowing the 
description of systems affected by nonlocal effects, such 
as long-range forces and memory effects. 

The study of astrophysical self-gravitating systems was 
one of the first applications of the theory [6j. This ap- 
proach has recently witnessed a renewed interest in the 
hope that it may provide a theoretical basis for the de- 
scription of the universal structure of dark-matter (DM) 
halos 0, H, H, H3|. Nonextensive statistical mechanics 
predicts their equilibrium states to be stellar polytropes 
(SP) [H, El, which have been claimed to fit DM halos as 
well as the usual Navarro, Frenk & White model ITol ] 
(hereafter NFW [Hj]). Moreover, SP have the distinct 
advantage of being analytically derived from nonexten- 
sive statistical mechanics, while most models describing 
DM halos are empirical fits to N-body simulations. In 
this context, nonextensive statistical mechanics appears 
as an attractive framework for providing a theoretical 
understanding of the structure of DM halos and self- 
gravitating systems in general. 

In this paper we compare, in a parameter indepen- 
dent way, the equilibrium configuration of astrophysical 
self-gravitating systems predicted by nonextensive sta- 
tistical mechanics to simulated DM halos. We clarify 
the issue of which central boundary conditions to use for 
comparing SP to simulated halos. On this basis, we es- 
tablish that simulated DM halos do not corroborate the 
predictions of nonextensive statistical mechanics. This 
result calls into question the direct applicability of nonex- 
tensive statistical mechanics to equilibrium astrophysical 
self-gravitating systems. 
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II. NONEXTENSIVE STATISTICAL 
MECHANICS AND STELLAR POLYTROPES 

Nonextensive statistical mechanics is based on the for- 
mulation of a generalized entropy 

S q = k B ^ (q&R), (1) 

which reduces in the limit q — > 1 to the Boltzmann-Gibbs 

entropy Sbg = —ks Y^Li Pi m Pi (P* 1S the probability of 
finding the system in the microstate i, W is the number of 
microstates i for a given macrostate, and ks is the Boltz- 
mann constant). When the index q ^ 1, the entropy of 
the system is nonextensive, i.e., it is not possible to sum 
the entropies S q (A) and S q (B) of two subsystems A and 
B, but the generalized entropy S q is conveniently pseudo- 
additive: S q (A+B) = S q (A)+S q (B)+(l-q)S q (A)S q (B). 
The third term accounts for coupling between the sub- 
systems, due to nonlocal effects. 

In this framework, astrophysical self-gravitating sys- 
tems at equilibrium are described by SP. This solution 
arises from extremizing the entropy S q under the con- 
straints of fixed total mass and energy as derived in [rlfllj] 
(see Ref. [l3[ about the choice of constraints), leading to 
a spherically symmetric isotropic system. The resulting 
distribution function depends on the parameter q and, 
via the identification n = 1/2 + — ^ [32], has the same 

form as the SP distribution function: /(e) oc e n ~ 3 / 2 for 
e > and /(e) — for e < 0, where n is the polytropic 
index and e the relative energy per unit mass |14| . 

In order to compare the predictions of nonexten- 
sive statistical mechanics to the results of N-body 
simulations, we compute radial profiles of fundamen- 
tal quantities of astrophysical self-gravitating systems: 
the matter density p(r), the logarithmic density slope 
dlog(p)/dlog(r), the integrated mass M(r), and the cir- 
cular velocity V c (r) (see Fig. 1). SP are solutions of the 
Lane-Emden equation (T^ |. depending on three free pa- 
rameters: the polytropic index n, the central density pa 
and the central velocity dispersion uq. The value of n 
(and q) is an intrinsic property of each system, reflecting 
its degree of nonextensivity, but has not been determined 
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FIG. 1: From top to bottom are presented the density, the 
logarithmic density slope, the integrated mass and the circular 
velocity. The profiles are scaled to r_2, the radius at which the 
slope equals — 2. Stellar polytropes predicted by nonextensive 
statistical mechanics are represented by the solid curves with 
the shades of green corresponding to different polytropic index 
n. Simulated dark-matter halos are represented by the NFW 
model (pink dashed curve), the 3D Sersic model (pink dotted 
curve; a = 0.17) and the Hernquist model (pink dash-dotted 
curve; a — 0.45). 



from first principles and is therefore a free parameter. SP 
with n > 3/2 are stable due to Antonov's stability crite- 
rion, i.e., d/(e)/de < [l4|. Other values of n are un- 
realistic: n = 3/2 corresponds to a distribution function 
independent of e, and n < 3/2 to a distribution function 
diverging at the escape energy, e = 0. Therefore, we com- 
pute numerical solutions of the Lane-Emden equation for 
n > 3/2, choosing values of n ranging from 2 to oo. We 
verified the consistency of our profiles using the analyti- 
cal solutions of the Lane-Emden equation, as well as the 
isothermal sphere which is the asymptotic solution when 
n — > oo [331 ] . 

Nonextensive statistical mechanics does not specify the 
boundary conditions for SP, which are needed to solve 
numerically the Lane-Emden equation. Simulated ha- 
los necessarily have a finite gravitational potential at the 
center, therefore we look for SP fulfilling this condition. 
For SP the density p is linked to the gravitational po- 
tential tp by p = est x ip n so we require solutions 
of the Lane-Emden equation with finite density at the 
center. Chandrasekhar proved that such solutions have 
dlog(p)/dlog(r) = at r = [15]. Therefore, the class of 
SP we can compare consistently to simulated halos have 
the natural boundary conditions at the center: p = po 
and dlog(p)/dlog(r) = 0. We present the corresponding- 
profiles in Fig. 1. In [UGH, the authors arbitrarily fix a 
non-zero logarithmic density slope near the center, cho- 
sen to match the NFW density profile (34|. This class 
of solutions must be treated with caution, as it does not 
ensure physical boundary conditions at the center. 



III. COMPARISON TO SIMULATED DM 
HALOS 

Cold dark-matter (CDM) halos resulting from N-body 
simulations provide an excellent way of testing the pre- 
dictions of nonextensive statistical mechanics. They are 
collisionless systems whose particles interact via long- 
range gravitational interaction only, allowing a direct 
comparison to theory. DM structures in N-body simula- 
tions evolve from primordial fluctuations via hierarchical 
clustering to form halos of universal shape, well described 
by the NFW model [13] ■ These halos formed through ac- 
cretion and mergers are sufficiently similar to those ob- 
tained by monolithic collapse, i.e., isolated halos, so that 
we do not take into account in this study the influence of 
formation history and cosmological environment on the 
final state of simulated halos [l7[ • 

Collisionless systems have very long collisional relax- 
ation timescales to reach the 'true' equilibrium, but ac- 
cess a stable quasi-equilibrium state faster (dynamical 
timescales are short compared to cosmological timescales, 
i.e., the age of the universe) through the processes of 
violent relaxation and phase- mixing [LM ]. The quasi- 
equilibria predicted by numerical cosmological simula- 
tions and by nonextensive statistical mechanics, respec- 
tively, may in either case be considered the most probable 
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state the system reaches. 

Many empirical models have been proposed as fitting 
the universal profile of CDM halos. We use in our study 
three of them to represent simulated halos. The NFW 
model is the more commonly used [l2[ . The Sersic model 
has been proposed recently as fitting, more accurately the 
inner part of high resolution halos [H, [l9[ . The Hernquist 
model provides a simple description of sclf-gravitating 
systems, for which all the profiles we use have analytical 
expressions [13] • As seen in Fig. 1, these profiles form 
a group with similar shapes, and we do not need to use 
actual simulations in our study. 

To compare simulated DM halos to SP, we select the 
radial range for which recent N-body simulations are ro- 
bustly resolved. Based on the high resolution simula- 
tions published in Navarro et al. [18| , this corresponds to 
0.1r_ 2 < r < 10r_ 2 , where r_ 2 is the radius at which 
the logarithmic density slope equals —2. To compare SP 
of all values of n with simulated halos, we need a scaling 
suitable for finite and infinite mass halos, circular velocity 
profiles with and without maximum, and density profiles 
with very different steepness: r_2 provides such a univer- 
sal scaling [35]. Moreover, it allows to take equally into 
account the inner and outer parts of the halo, respectively 
defined as having a slope shallower and steeper than —2 
fl2| . Finally, our study is independent of scaling param- 
eters: SP profiles are scaled depending on the choice of 
po and (To, but the shape depends only on n. There- 
fore, scaling to r_2 makes identical any profiles of same 
n but different po, 0o, keeping only information about the 
shape. All the profiles we present, as scaled to r_2, de- 
pend only on structural parameters: range of n for SP; 
a = 0.17 for the 3D Sersic profile [jj]; a = 0.45 for the 
Hernquist profile [20| : r s for the NFW profile [3q . In 
Fig. 1 we compare SP with simulated CDM halos, scal- 
ing to r_ 2 - 

In the outer parts, where r > r^2, SP have very differ- 
ent shapes depending on the polytropic index n. Profiles 
with n < 5 correspond to finite-mass halos, while n > 5 
polytropes have infinite mass, tending to the isothermal 
sphere when n — > oo. While finite-mass polytropes may 
appear more attractive from a physical point of view, 
they also provide worse fits to simulated DM halos, with 
an outer slope as steep as dlog(p)/dlog(r) < —5 at 
10r_2- For infinite-mass SP with larger values of the 
polytropic index, n = 16.5 (as found in Q) agrees with 
NFW and Sersic models, n = 10 with the Hernquist 
model. Therefore, some infinite-mass SP can provide a 
good description of the outer parts of simulated halos. 

In the inner parts, where r < r_2, SP share similar 
properties: they have a large core, the shape and the 
extent of which depend little on the value of n. This 
feature is in striking disagreement with the predictions 
of N-body simulations, which show steeper inner slopes. 
While this core structure itself has been advanced as an 
advantage of SP over the NFW profile solving the 
well-known cusp problem between simulations and ob- 
servations [2l| . it appears in Fig. 1 that this is not the 



case: the difference between SP and simulated halos in 
the density, density slope and mass profiles is as high as 
one order of magnitude at 0.1r_2- Therefore, for any 
polytropic index n, SP do not properly describe the inner 
parts of simulated halos. 



IV. DISCUSSION 

We have established here that the predictions of nonex- 
tensive statistical mechanics are not corroborated by sim- 
ulations of DM halos. Our results are based on readily 
observable quantities, i.e., the matter density itself, and 
profiles derived from it. These findings are in direct con- 
trast to previous works which found reasonable agree- 
ment because they either considered only the outer parts 
of simulated halos [8|, or used a non-zero density slope 
as initial condition near the center 0], or used values of 
n < 5 that lead to too steep density profiles in the outer 
parts [l(J. Support for our conclusions, however, comes 
from Barnes et al. (l6j who, even though fixing a non- 
zero density slope near the center to match the NFW 
density profile, found inconsistency between SP and DM 
halo velocity dispersion profiles. While we are cautious 
about the physical relevance of such profiles, it is of inter- 
est to emphasize that nor cored neither cuspy polytropes 
can describe properly simulated DM halos. 

These results imply a fundamental difference in the 
matter distribution of astrophysical self-gravitating sys- 
tems predicted by nonextensive statistical mechanics and 
by N-body simulations. Such a discrepancy raises three 
main questions. Is our comparison with simulated DM 
halos valid? Is nonextensive statistical mechanics the 
proper theory to describe collisionless long-range inter- 
action systems at equilibrium? And how does this study 
relate to observed DM halos? We briefly address these 
issues below. 

la. We study idealized systems which do not take into 
account complex effects in DM halos such as velocity 
anisotropy [22] or triaxiality 23]. However, the error we 
make by using isotropic and spherically averaged models 
to represent DM halos is small compared to the discrep- 
ancy we find between SP and simulated halos. Moreover, 
we do not address here general problems of statistical me- 
chanics of self-gravitating systems like, e.g., infinite mass 
(24| . so as to focus on the issues specific to nonextensive 
statistical mechanics. 

lb. N-body simulations depend on various parameters, 
such as the choice of the softening length and the number 
of particles, which can introduce numerical effects in the 
resulting halos. The softening of the gravitational force 
on small scale, introduced to avoid two-body interaction, 
creates a core at the center of the halo, approximately 
the size of the softening length. The number of particles 
fixes the maximum phase-space density resolved at red- 
shift zero, and if not large enough, a core appears due 
to the lack of particles in the center. Both effects lead 
to a shallower density profile in the center of halos if the 
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simulation is of low resolution [25(. Therefore, the dis- 
crepancy we observe between simulated profiles and SP 
is genuine and increases with the resolution of numerical 
simulations. 

lc. We compare theoretical statistical predictions to 
cosmological simulations, but the latter depend on ini- 
tial conditions and cosmological parameters. Initial con- 
ditions are fixed by the shape of the power spectrum of 
initial fluctuations P(k) oc fc" a , i.e., by the choice of the 
index n s . From CMB observations, we have an indication 
that n s — 0.958 ± 0.016 (26|, but tests on cosmological 
simulations proved that they depend very little on the 
choice of n s and on the cosmological model used [27| . 
Therefore, the universal profiles of simulated halos can 
be compared to statistical theories as a general predic- 
tion. 

Id. The choice of a scaling is necessary to compare 
DM halos models, but can bring a visualization bias. We 
checked if the use of r_2 leads us to overestimate the 
disagreement between nonextensive statistical mechanics 
and N-body simulations. However, adjusting the fit in 
the outer parts of the halo, e.g., scaling to the virial ra- 
dius, increases even more the discrepancy in the inner 
parts, while scaling profiles to fit well at smaller radii 
makes the discrepancy appear in the outer parts of the 
halo. 

2. Nonextensive statistical mechanics relies on three 
assumptions: a) the generalized entropy S q is the right 
form to describe long-range interaction systems; b) the 
system is at equilibrium, and its most probable state is 
given by the maximum entropy principle; c) S q is maxi- 
mized at fixed energy, leading to a system whose distribu- 
tion function depends only on the energy per unit mass, 
and has isotropic velocity dispersion. Though simulated 
halos show evidence of velocity anisotropy, the relation 
between density slope and velocity anisotropy [22j shows 
that it is isotropic in the center, where violent relaxation 
[HI (the process by which strong potential fluctuations 



efficiently drives the system towards equilibrium) is most 
effective. Therefore, the assumptions b) and c) hold true 
in the center. However, it is in the inner parts that 
the disagreement between SP and simulated halos is the 
strongest. Hence, we suggest that the assumption a), 
i.e., the choice of the generalized entropy S q motivated 
by nonextensive statistical mechanics, cannot be used to 
predict the equilibrium structures of simulated DM halos. 

3. Simulated DM halos are in good agreement with 
observations [29l |. except for the inner core found from 
spiral galaxy rotation curves |21|. Stellar polytropes, in- 
terestingly, have an inner core too, but the discrepancy 
we observe between SP and simulated DM halos is too 
large to explain the core of observed DM halos. Adding 
adiabatic contraction in simulations results in DM ha- 
los with steeper inner parts, which would not change our 
conclusions. 

In summary, we have established that nonextensive 
statistical mechanics 

0, S HH, a theory generalizing 
classical statistical mechanics and thermodynamics, does 
not describe the equilibrium state of astrophysical self- 
gravitating systems, as represented by DM halos formed 
in N-body simulations. 
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